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ABSTRACT 

We investigate consequences of the discovery that Fe II emission in quasars, one of 
the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. 
Eigenvector 1 correlates with the Eddington ratio L/L E ^ so that Fe ill H/3 increases 
as L/LEdd increases. We show that the "force multiplier", the ratio of gas opacity 
to electron scattering opacity, is ~ 10 3 — 10 4 in Fe II emitting gas. Such gas would 
be accelerated away from the central object if the radiation force is able to act on the 
entire cloud. As had previously been deduced, infall requires that the clouds have large 
column densities so that a substantial amount of shielded gas is present. The critical 
column density required for infall to occur depends on L/L^^, establishing a link 
between Eigenvector 1 and the Fe ll/H/3 ratio. We see predominantly the shielded face 
of the infalling clouds rather than the symmetric distribution of emitters that has been 
assumed. The Fe II spectrum emitted by the shielded face is in good agreement with 
observations thus solving several long-standing mysteries in quasar emission lines. 
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1. Introduction 



Principal component analysi s (PCA) of the emission -line spectra of quasars isolates several 
linearly independent components (|Boroson & Greenlll992|) . These are ascribed to distinct physical 
processes within these objects. Eigenvector 1, which is responsible for the largest part of the vari- 
ability within the samples, anti-correlates with L/L Ec id, t he ratio of the luminosity to the Eddington 
limit, while Eigenvector 2 represents the accretion rate (IBorosonll2002|) . It has been proposed that 
Eigenvector 1 represents the eff ects of radiation pressure controlling the column densities of the 
surviving line-emitting clouds (|Marconi et al.ll2008l : lNetzerll2009l: Dong et al.ll2009|) . This Letter is 
a specific, quantitative exploration of that idea. 

Fe II emission is one of the main spectroscopic signatures of Eigenvector 1. The physics of 
Fe II emission is complex wit h many levels contributing to the spectrum, which occurs i n a blended 
complex of muddled features dwills et alJl985ilVerner et al.ll 19991 lialdwin et alJ2004h . There has 
been only partial succ ess in reproducing th e observed Fe II emission in a realistic model of AGN 
emission-line regions ([Baldwin et al.ll2004l) . 



Gaskelll (|2009|) has recently discussed the considerable range of exis ting obse r vationa l ev 



idence indicating that QSO emission-line regions include infalling gas. iHu et al.l (|200 8a) and 
Hu et al.1 (|2008bl) recently found that Fe II emission is systematically redshifted relative to the 
QSO systemic velocities, and presented the reasons suggesting that this is because the Fe II emis- 
sion comes from infalling clouds. In a f orthcoming pape r (C. Hu et al. 2009, in preparation), we 
will further develop a model outlined by lHu et al.l (|2008b|) in which the Fe II emission comes from 
infalling clouds systematically viewed from their shielded faces. Such clouds must have a large 
enough column density to fall toward the black hole, despite L/L E dd ~ 10" 1 — 1. Here, we de- 
scribe a preliminary exploration of the basic spectroscopic properties of clouds viewed from their 
non-illuminated sides, which is an important question in its own right independent of the exact 
dynamical model. If indeed the Fe II emission does come from infalling clouds, our results have 
major consequences for the predicted emission-line spectrum and reveal one of the underlying 
drivers for Eigenvector 1 . 



2. Observed properties of Fe n emitters 



Hu et al.l (|2008a) and IHu et al.l (|2008bl) find that much of the optical Fe II emission comes 



from a redshifted intermediate-line region. The redshift of the Fe n-emitting region inversely cor- 
relates with the L/L E dd rat i° as indic ated by Eigenvector 1. The correlations are that, as the L/L E dd 
ratio increases, Fe ill H/3 increases (|Boroson 2002 . Boroson & Green 1992 ) and the Fe II redshift 
relative to the systemic rest frame decreases ( Hu et al.ll2008b|) . These relations are the subject of 
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this Letter. 



For simplicity, we consider two Fe II emission bands. The o ptical Fe II A4558 b and is the 
integrated Fe II emission between A 4434 and A4684, as defined by lFrancis et al.1 (I1992|) . UV Fe II 
emission forms a bl ended pseudo continuum with very few isolated Fe II lines. We follow the 
Baldwin et al.l (|2004l) definition of UV Fe II A 2445 as the integrated Fe II flux over the wavelength 
range AA2240- 2650. 

We consider the properties of a typical Fe n-emitting cloud. The observed spectrum is actually 
produced by a mix o f clouds with a broad range of densities and distances from the central object 
([Baldwin et al.lll995l) . However, in this preliminary investigation, we consider only a single cloud. 



We use a combination of line-continuum r everberation studie s and theoretical predictions to 
determine the parameters of this typical cloud. iMaoz et al.1 (|1993[) found that the UV Fe II lines 
in NGC 5548 respond to the driving continuum on a timescale similar to that of C IV A 1549 and 
somewhat shorter than that of H/ 3, but the optical Fe II lines in t his object are too w eak for use- 
ful reverberation measurements (lVestergaard & Peterson! |2005|) . iKuehn et al.1 (|2008[) studied the 
reverberation behavior of the optical Fe II lines in Akn 120, for which there are no UV rever- 
beration measurements but which has strong optical Fe II lines. They found that in Akn 120 the 
optical Fe II emission clearly does not originate in the same region as H/3, and that there was some 
evidence of a reverberation response time of 300 days which implies an origin in a region sev- 
eral times further away from the central object than H/3 . A 300 day light-travel time corresponds 
to a source-cloud separation of ~ 8 x 10 17 cm. For a general AGN spectral energy distribution 
(SED), the quoted luminosity of Akn 120 (10 45 erg s _1 ), and this separation, the flux of hydrogen- 
ionizing photons, 4>(H), at the derived position of the Fe II e mission is in the nei ghborhood of 



<3>(H) pa 10 19 cm 2 s 1 . From the results shown in Figure 7 of IBaldwin et al.l (I2004|) . we estimate 



the gas density to be «(H) ~ 10 10 cm 3 and the turbulence to be u tur t, ~ 10 2 km s 1 . The "standard 
cloud" considered in the remainder of this Letter has these quantities and solar abundances. 



3. Calculations 



All calc ulations were done with version C10 of the spectral simulation code Cloudy, last 
described by IFerland et al.l (| 19981) . We concentrate on the emission from the shielded face of the 
cloud. Cloudy has long pre dicted the fraction of the emission that is beamed toward the central 
object (|Ferland et al.l Il992|) . but the outward emission, which is often quite faint for the large 
column densities we consider below, has not been an emphasis. The inward versus outward Fe II 
emission, although calculated, was not reported. We have improved the adaptive logic used to 
define the spatial grid to better resolve the emission from the shielded face and added several 
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output options to report these predictions . We now exp l icitly predict the inward, outward, and total 
Fe II emission in the form presented by lBaldwin et al.l ([2004). 



We consider the ionizing flux versus gas density plane presented by iBaldwin et al.l (|2004|) 



For an isotropic radiation field, the ionizing flux ^>(H) 



so this variable can be thought of as 



a proxy for the continuum source-cloud separation. All possible AGN clouds lie somewhere on 
this plane. Regions with the same ionization parameter form parallel lines with a slope of unity. 
Little observable emission comes from the large flux, small de nsity quadrant whic h corresponds to 
highly ionized gas with temperatures near the Compton limit (jKorista et al.ll 19971) . Clouds in the 
small flux, large density quadrant have very low ionization, tend to be cool, and produce strong 
Fe II emission. 



3.1. The force multiplier 



The fact that the Fe n-emitting clouds, with L/LEdd ~ 10 , are falling into the center has 
major consequences for the cloud column density. The Eddington limit is the luminosity at which 
the force due to electron scattering balances the gravitational acceleration for optically thin gas, 



-•Edd 



AnGMnipC 



(1) 



Here Gt is the Thomson cross section. The ratio L/LEdd is a dimensionless measurement of the 
specific luminosity of the AGN. 

This acceleration limit is appropriate for very highly ionized gas. For gas with a higher total 
cross section a, the outward force is greater by a force multiplier g/gt- In general, G ^> Gt for 
low-ionization gas due to the added opacity cau sed by line and photoelectric absorption (|Castor 
19741 : iKippenhahn et al.ll 19741 : ICastor et al.lll975i) . The force multiplier over the range of densities 
and ionizing fluxes i s shown in Figure [H which also shows, as a shaded ellipse, the region that 
Baldwin et al.l (|2004) found could produce acceptable Fe II emission. 



Figure \T\ shows that typical Fe il-emitting clouds should be strongly accelerated away from 
the centr a l objec t unless the Eddington ratio is below L/L E dd ~ 10 3 — 10 4 . The objects in the 
Hu et al.l (|2008bh sample have L/LEdd ~ 10 _1 so the gas would be expected to be accelerated 
outward. How can it be in a state of infall? 
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Fig. 1. — Force multiplier, evaluated at the illuminated face of the cloud, over the photon flux- 
density plane. The s haded ellipse is the region capable of reproducing Fe II observed emission 
JBaldwin et al]l2004h . 
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3.2. The need for large column densities 



The outward radiation pressure corresponding to the force multiplier shown in Figure Q] will 
act only on the illuminated face of a large column density cloud. This is because the peak of the 
SED emitted by the AGN occurs at ionizing energies. The photons most capable of pushing the 
gas away ionize the gas and are absorbed. 

An H + -H° ionization front occurs when most ionizing photons have been absorbed. Most of 
the outward momentum acts on the H + layer. Relatively little energy is in high-energy photons 
which penetrate into neutral regions. The result is that most of the radiative acceleration occurs 
in the H + layer which then pushes on deeper neutral and molecular regions. The presence of 
shielded regions allows inflow t o occur even at super- Eddi ngton luminos ities, as has been discussed 
previously (e.g.. IShavivll 1 998|) . iMarconi et al.l (I2008|) and lNetzerl (|2Q09r> have previously discussed 
radiation pressure from ionizing photons acting on large column-density AGN clouds, in terms of 
the effect on the deduced black hole masses and the L/L^^ versus Mbh relation. 

We quantify this with a calculation shown in Figured using the standard conditions described 
in Section [2] The gas kinetic temperature is shown as a function of column density from the 
illuminated face of the cloud. A sharp drop in temperature occurs at the H + -H° ionization front 
and an H + layer with a column density of N(H) ~ 5 x 10 21 cm~ 2 exists on the face of the cloud. 
The vast majority of the cloud consists of warm atomic gas where the Fe II lines form. 

The dashed line shows the computed radiative acceleration. Within the H + layer the force 
increases nearly linearly with increasing column density as more of the incident radiation field is 
absorbed. The force has nearly reached its asymptotic value when the H + -H° ionization front is 
reached. This ionized layer then pushes against the much larger column density of cooler atomic 
gas. The atomic layer adds mass to the cloud without contributing to its outward acceleration. 

In practice, the clouds may be subject to Rayleigh-Taylor instability ( Mathews1ll982 , Mathews 



1986|) . While the radiation force will suppress instability at the front surface of the cloud, it will 
enhance the instability at the ionization front. Within the ionized region, the struc ture may be 
subject to additional radiation-driven instabilities (fMestel et al.lll976l : IWilliamsl 120001) . as well as 
significant structural perturbations due to the varying continuum. These and other sources may 
support a level of turbulence similar to that which we model. We will not address the stability of 
the clouds we model further. The model does require that the cloud integrity must be maintained 
by other means, possibly magnetic, or, more likely, that the cloud is an evolving dynamic entity 
that survives long enough to be accelerated by the gravitational forces in the region. 

Figure [2] establishes the minimum column density where infall is possible. The most useful 
way to consider the results in Figure |2] is by reference to the acceleration that would occur for 
pure electron scattering opacity in fully ionized gas since that is the opacity used in the definition 
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Fig. 2. — Gas kinetic temperature, the computed force for the radiative acceleration, and the force 
for optically thin electron scattering, are shown as a function of column density from the illumi- 
nated face. Ionizing radiation enters from the left. The two lines showing the integrated force 
intercept at a column density of N(H) ~ 1.2 x 10 24 cm~ 2 , the point where the cloud would have 
neutral buoyancy for an AGN at the Eddington Limit. 
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of the Eddington limit. The dash-dotted line shows the equivalent radiative driving for the case 
where only electron scattering opacity acting on the incident radiation field is considered. The 
actual acceleration, the dashed line, is considerably larger than the electron scattering acceleration 
across the H+ region, as expected from Figured] The acceleration increases very slowly after the 
H + - H° ionization front, while, with our definitions, the reference electron scattering acceleration 
increases linearly with column density Af(H). The two lines cross at a column density of N(H) ~ 
1.2 x 10 24 cm~ 2 , the point where the cloud would be neutrally buoyant for L/L,Edd = 1- 

The minimum column density for which the inward gravitational force is stronger than the 
total outward radiation pressure on a cloud, iV(H)infall> ma Y be estimated as 

^(H) inf aii > ^ L5 X ^VtMck^-cm- 2 , (2) 

&T L-Edd L Edd 

where /thick ~ 1 is the fraction of the incident radiation to which the cloud is optically thick. This 
establishes a relation between specific luminosity and the minimal column density of infalling 
clouds. In objects of higher specific luminosity, which tend to have stronger Fe II emission, only 
high column density clouds will be able to fall inward. 

Figure [3] shows how the total and, separately, the outward emission varies as the cloud col- 
umn density is increased. The Fe II lines begin to form when clouds have column densities 
N(H) > 10 cm~ 2 , the point where the H + -H° ionization front occurs. This is also the column 
density where H/3 becom es inwardly beamed. La is nearly fully emitted in the inward direction 



(Ferland_&NetzeJlL979) . 



3.3. Isotropy of optical Fe n emission 



The upper panel of FigureHJshows the total and outward Fe II emission for our standard cloud. 
The spectrum is complex and it is hard to distinguish between the total and outward components 
in certain parts. The lower panel shows the ratio of the outward to the total emission. This shows 
the general trend for the inward fraction to increase with decreasing wavelength, mainly due to the 
larger optical depths of UV Fe II transitions. 

Figu re [5] compares the outward e mission from our standard mod el with the Fe II template pre- 
sented by IVestergaard & Wilkesl (1200 lb and lBoroson & Greenl (|1992l) . The agreemen t is good. We 
notice in particular that the predicted optical / UV ratio, which in the simulations by lBaldwin et al. 



(120041) is always too small when both inward and outward Fe II emission is considered, is actually 
larger than observed (i.e., the UV Fe II strength is underpredicted). This allows for a component 
of UV Fe II emis sion to form in the broad-line region (BLR), as suggested by reverberation results 
of NGC 5548 dlVlaoz et all 19931) . 
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Fig. 3. — Predicted total (solid line) and outward (dashed line) emission relative to the total H/3 
is shown as a function of column density. H/3 is isotropically emitted for small column densities 
while the optical Fe II band remains nearly isotropic for all column densities. The UV Fe II band 
is predominantly inwardly beamed due to the larger optical depth in these lines. 
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Fig. 4. — Upper panel compares the total (solid) and outward (dotted) Fe II emission of our stan- 
dard cloud. The lower panel shows the ratio of the outward to total emission. 
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This is a major advance in the comparison between theory and observations for the Fe II 
emission. It is important to keep in mind the fact that the deduced properties of emission-line 
clouds have, until now, come from comparing the total emission, not the outward emission, with 
observations. Future papers will further investigate the implications of asymmetric anisotropic 
emission. 



4. Discussion 



Equation © suggests a linear relationship between the Eddington ratio and the column density 
of infalling clouds. Figure[3]shows that the Fe ll/H/3 ratio increases with increasing column density 
for N(R) ~ 10 22 — 10 23 cm~ 2 . This implies that there will be a relation between L/Lgdd an d 
Fe ll/H/3. This explains one of the strongest syste matic correlations captured by Eigenvector 1. As 
was also recently suggested by lDong et al.1 (|2009), the cloud column density is the physical variant 
which drives the spectroscopic variations. 

The predicted Fe ll/H/3 reaches an asymptotic value of about 3 for large Af(H) or L/L^d- This 
is close to the observed value for luminous quasars ([Francis et al.lll99lf). Models which considered 
the full emission from the cloud do not reproduce this ratio ([Baldwin et a l. 2004). 



La is strongly inwardly be a med, while C IV i s nearly isotropic, as shown by lFerland & Netzei 



(119791) and lFerland et all (1 19921) . iHu et all (I2008bl) find that all of the emission in the optical Fe II 
band originates in the infalling clouds. Other lines are likely to form in both these and other 
regions, including the accretion disk, virialized broad-lined clouds and outflowing winds. High- 
resolution UV spectra of objects in which the infalling clouds have the largest velocity shift are 
needed to test this. 

The infalling clouds must have a substantial coverin g factor. The o bserved H/3 equivalent 
widths of the broad and infalling components are similar (|Hu et al.ll2008al) . The infalling clouds 
direct much of their H/3 emission away from the observer. They must have a coverin g factor 
approaching Q./Atz ~ 0.5, consistent with an origin in the molecular torus (|Hu et al.ll2008b|) . 



The net accretion rate due to these infalling clouds must be 



A^Fell > -r^r l m p {N H /r)vmm 
4-7C 



a 

17— ri 8 Af H ,24V8M yr 
in 



(3) 



where r = 10 18 riscm, Nu = 10 24 /VH.24cm~ 2 , and v = lOOOvskms -1 . In this expression, the factor 
Nn/r is the mean density of the Fe Il-emitting material multiplied by the minimal line-of- sight 
filling factor if the clouds are evenly distributed. This may be compared to the accretion rate of 
1.8L45(t]/0.01) _1 M yr _1 required to power an AGN of luminosity L = lO^Z^ergs -1 with an 
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accretion efficiency of 77, where the scaling parameters are all ~ 1 for Akn 120. The mass infall 
rate above that required by accretion is likely to be balanced by other outflows. Nevertheless, this 
suggests that the infalling Fe n-emitting material constitutes a major element in the mass budget 
of high Eddington-ratio active nuclei. 

We have investigated the effects of systematically viewing a low-ionization cloud from its 
shielded face. We focus on Fe II and H/3 because of their importance in PCA analysis. Meaningful 
predictions about the full spectrum will require a much broader exploration of, and probably a 
weighted integration over, the ionizing flux versus gas density plane. A future paper will carry out 
that broader exploration and address predicted strengths and profiles of such lines, and also the 
open question of whether this infalling component should be seen in absorption particularly at UV 
or X-ray wavelengths. 



Conclusions 



• Low-ionization Fe il-emitting clouds have force multipliers (the ratio of total to electron 
scattering opacities) that are large if the clouds are optically thin. Such clouds will be accel- 
erated away from the central black hole if it has an Eddington ratio as large as those found 
in quasars, L/L Edd ~ 0.1. 

• If the Fe il-e mitting clouds are infalling, as is suggested by their systemically positive ve- 
locity offsets (Eu et al. 2008b), they must have column densities significantly larger than 
previously thought if the inward pull of gravity is to offset the outward force of radiation 
pressure. 

• There is a simple relationship between L/L Edc j and the minimum cloud column density re- 
quired for infall. Fe ll/H/3 also depends on cloud column density. Cloud column density is 
the underlying driver that c ouples Eigenvecto r 1, L/L^^, and spectroscopic variations, as 
was previously proposed by Dong et al.1 (|2009h . 



• Previous simulations of the emitting gas have assumed a symmetric distribution of clouds 
so that, on average, we see the same number of clouds from their illuminated as from their 
shielded faces. We have investigated here the result of seeing an asymmetric distribution 
of emitters, so that we mainly observe Fe II emi ssion from the shielded face of near-side 
infalling clouds, as is suggested by observations (|Hu et al.ll2008bl iGaske ll 2009). 



• Fe II UV emission is emitted less isotropically than optical Fe II lines. Previous work, which 
assumed a symmetric distribution of emitters, could not explain the larger-than-predicted 
ratio of optical to UV Fe II emission or the Fe ll/H/3 ratio. The predicted emission ratios 
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from the shielded face are in good agreement with observations. This is a major advance in 
understanding the nature of Fe II emission in AGN. 

It is likely that there is a distribution of cloud column densities. If so, clouds with column 
densities smaller than that given in Equation © will be radiatively accelerated outward while 
those with larger column densities will fall in, with a column density cutoff dependi ng on the 
AGN ' s luminosity. Th is would be the mechanism that accounts for the observed (|Boroson 
2002l : lDong et al.ll2009|) coupling between L/L^dd and the column density and Fe ll/H/3 ratio. 
Both outflow and infall are observed in AGN. Could this also be a natural consequence of a 
range of N(R) in the surviving clouds? 
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Fig. 5. — Red smoother line is the outward Fe II emission from the standard model broadened to a 
FWHM of 900 km s while the black spiky line is the observed Fe II template. Intermediate-line 
region clouds account for the optical emission while allowing part of the UV emission to come 
from the inner broad-line region. 



